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I. INTRODUCTION 

d\ '. 

The 0(3) nonlinear sigma model has long been the subject of intense research due to its theoretical and phenomeno- 
logical basis. This theory describes classical (anti) ferromagnetic spin systems at their critical points in Euclidean 
space, while in the Minkowski one it delineates the long wavelength limit of quantum antifcrromagncts. The model 
O ' exhibits solitons, Hopf instantons and novel spin and statistics in 2+1 space-time dimensions with inclusion of the 
, Chern-Simons term. 

The soliton solutions of the model exhibit scale invariance which poses difficulty in the particle interpretation on 
^■f-^ , quantization. A popular means of breaking this scale invariance is to gauge a U(l) subgroup of the 0(3) symmetry 
^-H ' of the model by coupling the sigma model fields with a gauge field through the corresponding U(l) current^! This 
[ class of gauged 0(3) sigma models in three dimensions have been studied over a long time^'^'^i^'^'^'i. Initially the 
. gauge field dynamics was assumed to be dictated by the Maxwell terml. Later the extension of the model with the 
^ ' Chern - Simons coupling was investigated^. A particular form of self - interaction was required to be included in 
these models in order to saturate the Bogomol'nyi bounds^. The form of the assumed self - interaction potential is 
of crucial importance. The minima of the potential determine the vaccum structure of the theory. The solutions 
change remarkably when the vaccum structure exhibits spontaneous breaking of the symmetry of the gauge group. 
I Thus it was demonstrated that the observed degeneracy of the solutions is lifted when potentials with symmetry 
On , breaking minima were incorporated^i^. The studies of the gauged 0(3) sigma model is important due to their intrinsic 
On > interest and also due to the fact that the soliton solutions of the gauged 0(3) Chern-Simons model may be relevant 
in planar condensed matter systems^iiSiiii. Recently gauged nonlinear sigma model was considered in order to obtain 
. self-dual cosmic string solution s '^^i^^ . This explains the continuing interest in such models in the literature-i^ii^ii^i^i^ii. 
O i' A particular aspect of the gauged 0(3) sigma models where the gauge field dynamics is governed by the Maxwell 
, term can be identified by comparing the results of^ and^. InX the vaccum is symmetric and the n — 1 soliton 
' solution does not exist, n being the topological charge. Here, solutions exists for n = 2 onwards. Moreover, these 
^ , soliton solutions have arbitrary magnetic flux. When we achieve symmetrybreaking vaccum by chosing the potential 
appropriatel}*^ soliton solutions are obtained for n ~ l.^^These solutions have quntized magnetic flux and qualify as 
magnetic vortices. It will be interesting to follow the solutions from the symmetrybreaking to the symmetric phase. 
^ . This is the motivation of the present paper. 
■ " ' We will consider a generalisation of the models oft and^ with an adjustable real parameter v in the expression of 
the self - interaction potential which interpolates between the symmetric and the symmetrybreaking vaccua. This 
will in particular allow us to investigate the soliton solutions in the entire regime of the symmetrybreaking vacuum 
structures and also to follow the collapse of the n — \ soliton as we move from the assymmetric to the symmetric 
phase. The solitons of the model are obtained as the solutions of the self - dual equation obtained by saturating the 
Bogomoln'yi bounds. Unfortunately, these equations fall outside the Liouville class even after assuming a rotationally 
symmetric ansatz. Thus exact analytical solutions are not obtainable and numerical methods are to be invoked. 

The organisation of the paper is as follows. In the following section we present a brief review of the 0(3) nonlinear 
sigma model. This will be helpful in presenting our work in the proper context. In section 3 our model is introduced. 
General topological classifications of the soliton solutions of the model has been discussed here. In section 4 the 
saturation of the self-dual limits has been examined and the Bogomol'nyi equations have been written down. Also 
the analytical form of the Bogomol'nyi equations has been worked out assuming a rotationally symmetyric ansatz. 
These equuations, even in the rotationally symmetric scenario, are not exactly integrable. A numerical solution has 
been performed to understand the details of the solution. A fourth order Runge - Kutta algorithm is adopted with 
provision of tuning the potential appropriately. In section 5 the numerical method and some results are presented. 
We conclude in section 6. 
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II. 0(3) NONLINEAR SIGMA MODELS 

It will be useful to start with a brief review of the nonlinear 0(3) sigma modeli^. The lagrangian of the model is 
given by, 



C = • (1) 



Here ^ is a triplet of scalar fields constituting a vector in the internal space with unit norm 

0a = na-</-, (a = 1,2,3) (2) 

(p ■ 4> = 4>a4>a = 1 (3) 

The vectors ria constitute a basis of unit orthogonal vectors in the internal space. We work in the Minkowskian space 
- time with the metric tensor diagonal, Qf^u = (1,-1,-1). 

The finite energy solutions of the model ([T]) satisfies the boundary condition 

Yimcj,'^ = 0^0) (4) 

at physiacal infinity. The condition ^ corresponds to one point compactification of the physical infinity. The physical 
space i?2 becomes topologically equivalent to 5*2 due to this compactification. The static finite energy solutions of the 
model are then maps from this sphere to the internal sphere. Such solutions are classified by the homotopyi^ 

n2(52) = z (5) 

We can construct a current 

1 

which is conserved irrespective of the equation of motion. The corresponding charge 



K^, = — e,,.A0 • {d''4> X aV) (6) 



1 

8^ 



d2xey<?!>-(9Vx'^'^) (7) 



gives the winding number of the mapping ([5]^ 



III. OUR MODEL - TOPOLOGICAL CLASSIFICATION OF THE SOLITON SOLUTIONS 

In the class of gauged models of our interest here a U(l) subgroup of the rotation symmetry of the model ([T]) is 
gauged. We chose this to be the S0(2) [U(l)] subgroup of rotations about the 3 - axis in the internal space. The 
Lagrangian of our model is given by 



C = -D^cj) ■ D^^cP - -F^.F^'^ + U{4>) (8) 
Df^iCj) is the covariant derivative given by 

Df^cf} = d^,cj) + A^n3 X (9) 

The S0(2) (U(l)) subgroup is gauged by the vector potential whose dynamics is dictated by the Maxwell term. 
Here F^j, are the electromagnetic field tensor, 

F^. = a^A, - d^Af, (10) 
U{<p) is the self - interaction potential required for saturating the self - dual limits. We chose 

u{4>) = -liv-<p3f (11) 
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where z; is a real parameter. Substituting v = we get back the model whereas v = 1 gives the model o^. 
We observe that the minima of the potential arise when , 



which is equivalent to the condition 



(12) 



02 + 02 ^ 1 - ^;2 (13) 

on account of the constraint ^ . The values of v must be restricted to 

\v\ < 1 (14) 

The condition (|13p denotes a latitudinal circle (i.e. circle with fixed latitude) on the unit sphere in the internal space. 
By varying v from -1 to +1 we span the sphere from the south pole to the north pole. It is clear that the finite energy 
solutions of the model must satisfy (I13p at physical infinity. For v 1 this boundary condition corresponds to the 
spontaneous breaking of the symmetry of the gauge group and in the limit 

\v\ ^ 1 (15) 

the asymmetric phase changes to the symmetric phase. We call the potential (jlip interpolating in this sense. In the 
asymmetric phase the soliton solutions are classified according to the homotopy 

ni(5i) = Z (16) 

instead of ([5]). In the symmetric phase, however, this new topology disappears and the solitons are classified according 
to (O as in the usual sigma model ([T]). A remarkable fallout of this change of topology is the disappearence of the 
soliton with unit charge. The fundamental solitonic mode n = I (n being the vorticity) ceases to exist in the symmetric 
phase. The modes corresponding to n = 2 onwards still persist but the magnetic flux associated with them ceases to 
remain quantised. 

In the asymmetric phase the vorticity is the winding number i.e. the number of times by which the infinite physical 
circle winds over the latitudinal circle Associated with this is a uniqe mapping of the internal sphere where the 
degree of mapping is usually fractional. By inspection we construct a current 

= ^eM-A[0 ■ D'^cj, X D^cl) - F-^{v - ^3)] (17) 
ovr 

generalising the topological current The current p7p is manifestly gauge invariant and differs from ([5]) by the 
curl of a vector field. The conservation principle 

d,.K^ = (18) 
thus automatically follows from the conservation of ([6]). The corresponding conserved charge is 



Using (jlTp and ((T9)l we can write 



T = y d?xKo (19) 
{v - h)A0rd9 (20) 



47r 



boundary 



where rfi are polar coordinates in the physical space and Ag = eg ■ A. Using the boundary condition ([T^ we find 
that T is equal to the degree of the mapping of the internal sphere. Note that this situation is different froniS. where 
the topological charge usually differs from the degree of the mapping. In this context it is interesting to observe that 
the current (|17p is not unique because we can always add an arbitrary multiple of 



with it without affecting its conservation. We chose ([TT]) because it generates proper topological charge. 
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IV. SELF DUAL EQUATIONS IN THE ROTATIONALLY SYMMETRIC ANSATZ 



In the previous section we have discussed the general topological classification of the solutions of the equations of 
motion following from ([5]) . In the present section we will discuss the solution of the equations of motion. The Euler - 
Lagrange equations of the system ^ is derived subject to the constraint ^ by the Lagrange multiplier technique 



where 



Using (1^ we get 



= [D,{D''ct,)-cf>]cP + n3{v~cl>3) 



-ng • and = x 



i^^J" = -(w-(/)3)(n3 X 
From (|22|) we find, for static configurations 

From the last equation it is evident that we can chose 

^0 = 



(21) 
(22) 



(23) 
(24) 
(25) 
(26) 



As a consequence we find that the excitations of the model are electrically neutral. 

The equations (PT|) and are second order differential equations in time. As is well known, first order equations 
which are the solutions of the equations of motion can be derived by minimizing the energy functional in the static 
limit. Keeping this goal in mind we now construct the energy functional from the symmetric energy - momentum 
tensor following from ([8|). The energy 



E 



1 



For static configuration and the choice — 0^ E becomes 

1 



E = 



d'x[{D,4,) ■ {D,cj)) + Fl^ + {v ~ 



Several observations about the finite energy solutions can be made at this stage from 

= 01 + i(l>2 

we get 

• =\{^^ + iA^)^f + {d^(l)3f 

The boundary condition (fT^l) dictates that 

4> « (1 -i;2)^e™e 

at infinity. From (|28p we observe that for finite energy configurations we require 

A = ee- 



By defining 



(27) 

(28) 

(29) 
(30) 
(31) 
(32) 



on the boundary. This scenario is exactly identical with the observations of^ and leads to the quantisation of the 
magnetic flux 



$ = / Bd^x = 



AgrdO = 27 



(33) 



boujidary 
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The basic mechanism leading to this quantisation remains operative so far as v is less than 1. At f = 1, however, the 
gauge field A becomes arbitrary on the boundary except for the requirement that the magnetic field B should vanish 
on the boundary. Remember that not all the vortices present in the broaken phase survives this demand. Specifically, 
the n = 1 vortex becomes inadmissible. 

Now the search for the self - dual conditions proceed in the usual way. We rearrange the energy functional as 

Equation ([M)) gives the Bogomol'nyi conditions 

Di(f)±eij(j)X Dj(jy = (35) 
Fi2 T (f - 03) - (36) 

which minimize the energy functional in a particular topological sector, the upper sign corresponds to +ve and the 
lower sign corresponds to -ve value of the topological charge. 

We will now turn towards the analysis of the self - dual equations using the rotationally symmetric ansata^^ 

(j)i{r,9) = sin g{r) cos n9 
4'2{t,9) — sin g{r) sin nO 
= cos g{r) 

AM = -ee^ (37) 
r 

From (jl2p we observe that we require the boundary condition 

g{r) cos^^ V as r — > oo (38) 

and equation ([32]) dictates that 

a(r) — 1 as r — > c» (39) 

Remember that equation (|32p was obtained so as the solutions have finite energy. Again for the fields to be well 
defined at the origin we require 

g{r) ^ or TT and a(r) -^0 as r — > (40) 
Substituting the Ansatz ([57|) into ([55)1 and ([55]) we find that 

5'(r)=±^ ^sing, (41) 

r 

a'(r) = ±-(i; - cos 5) (42) 
n 

where the upper sign holds for +ve T and the lower sign corresponds to -ve T. Equations (|4ip and ((42|) are not exactly 
integrable. In the following section we will discuss the numerical solution of the boundary value problem defined by 
(|iT|) and g21) with §^ to 

Using the Ansatz p7|) we can explicitly compute the topological charge T by performing the integration in (fT9|) .The 
result is 

n 1 

T = --[cosg(c») - cos.g(0)] - ~ cosg(c»)] (43) 

The second term of vanishes due to the boundary condition Also, when g(0) = 0, 

T=^il-v) (44) 

and, when g(0) = tt, 

T^-^{l + v) (45) 
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It is evident that T is in general fractional. Due to ([20)l it is equal to the degree of mapping of the internal sphere. 
This can also be checked explicitly. 

From the above analysis we find that g(0) — corresponds to +ve T and g(0) = tt corresponds to -ve T. We 
shall restrict our attention on negetive T which will be useful for comparision of results with those available in the 
literature. The boundary value problem of interest is then 

g'{r)^ ^ ^smg (46) 

r 

a'{r) = --{v - cosg) (47) 
n 

with 

g{0) = 7r,a(0) = 

(ji(cx)) = cos^"^ w, a(oo) = (48) 

In addition we require a'(r) ^ as r — > oo. This condition follows from (|¥71) and (HS)) and should be considered 
as a consistency condition to be satisfied by their soloutions. 



V. NUMERICAL SOLUTION 



The simultaneous equations and (|Tf|) subject to the boundary conditions ([35]) are not amenable to exact 
solution. They can however be integrated numerically. We have already mentioned the quenching of the n ~ 1 
solution in the limit v 1. This is connected with the transition from the symmetry breaking to the symmetric 
phase. The numerical solution is thus interesting because it will enable us to see how the solutions change as we 
follow them from the deep assymetric phase {v = 0) to the symmetric phase {v ~ 1). In the following we provide the 
results of numerical solution to highlight these issues. 

Let us note some details of the numerical method. A fourth order Runge - Kutta method was employed. The point 
r = is a regular singular point of the equation. So it was not possible to start the code from r = 0. Instead, we 
start it from a small value of r. The behaviour of the functions near r = can be easily derived from and (|Tf)) 

g{r) « TT + Ar" (49) 



a(r)«-^(l + «) (50) 

Here A is an arbitrary constant which fixes the values of g and a at infinity. In the symmetrybreaking phase the 
numerical solution depends sensitively on the value of ^.^^There is a critical value oi A, A = j^cnt j-^j. .y^rhich the 
boundary conditions are satisfied. If the value of A is larger than A'^"* the conditions at infinity are overshooted, 
whereas, if the value is smaller than the critical value g(r) vanishes asymptotically after reaching a maximum. The 
situation is comparable with similar findings elsewhere^. The values of g and a were calulated at a small value of r 
using and ([50)1 . The parameter A was tuned to match boundary conditions at the other end. Interestingly, this 
matching is not obtainable when n = 1 and v = 1. This is consistent with the quenching of the n = 1 mode in the 
symmetric vacuum situation. 

After the brief discussion of the numerical technique we will present a summary of the results. As may be recalled, 
the purpose of the paper is to study the solutions throughout the asymmetric phase with an eye to the disappearence 
of the n = 1 mode. Accordingly profiles of g and a will be given for n = 1, for different values of v. In figures 1 and 
2 these profiles are shown for t; = 0, .2, .4, .6, .8. The corresponding magnetic field distributions are given in figure 3. 
Another interesting issue is the change of the matter and the gauge profiles with the topological charge. In figures 4 
and 5 this is demonstrated for different n values for a constant v. 



FIG. 1: The matter profile. The function g[r) is plotted against r for n = 1 and different v values, indicated on the right hand 
top corner. 




2 4 6 8 10 



FIG. 2: The gauge field profile. The function a(r) is plotted against r for n — 1 and different v values. 




VI. CONCLUSION 

The 0(3) sigma model in (2+1) dimensional space - time with its U(l) subgroup gauged was mootedl as a possible 
mechanism to break the scale - invariance of the soutions of the original 3 - dimensional 0(3) sigma model. The 
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model finds possible applications in such diverse areas such as planar condensed matter physics^ii^iii, gravitating 
cosmic stringsiiiiS and as such is being continuously explored in the literatureii^i^ii^i^i^iiiii. An interesting aspect of 
the gauged 0(3) sigma models is the qualitative change of the soliton modes in the symmetric and symmetry breaking 
vacuum scenario, as can be appreciated by a comparision of solutions given ii)^. In this paper we have considered a 
gauged 0(3) sigma model with the gauge field dynamics determined by the Maxwell term as in-'^'^. An interpolating 
potential was included to invesigate the solutions in the entire symmetrybreaking regime . This potential depends 
on a free parameter, the variation of which effects transition from the asymmetric to symmetric phase. We have 
discussed the transition of the associated topology of the soliton solutions. The Bogomol'nyi bound was saturated to 
give the self - dual solutions of the equation of motion. The self - dual equations are, however, not exactly solvable. 
They were studied numerically to trace out the solutions in the entire asymmetric phase with particular emphasis on 
the n = 1 mode. Our analysis may be interesting from the point of view of applications, particularly in condensed 
matter physics. 
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This is different from the minimal coupling via the topological current discussed previously^". 

The disappearence of the n — 1 soliton has been shown to follow from general analytical method ii>^. 

This should be contrasted with the symmetric vacuum solution where A may have arbitrary values. 



